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Abstract 

We discuss how a dominant tetraquark component of the lightest scalar mesons may emerge in 
AdS/QCD gravity duals. In particular, we show that the exceptionally strong binding required to 
<^ ■ render the tetraquark ground state lighter than the lowest-lying scalar quark-antiquark nonet can 

(n: 

be holographically encoded into bulk-mass corrections for the tetraquark's dual mode. The latter 
O 

are argued to originate from the anomalous dimension of the corresponding four-quark interpolator. 
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To provide a concrete example, we implement this mechanism into the dilaton soft-wall dual for 



holographic QCD. Preventing the lowest-lying dual mode from collapsing into the AdS boundary 

JL , then establishes a rather generic lower bound on the tetraquark mass (which may be overcome in 

the presence of additional background fields). We further demonstrate that the higher tetraquark 



excitations can become heavier than their quark-antiquark counterparts and are thus likely to 
dissolve into the multiparticle continuum. 
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^ . PACS numbers: 11.25.Tq, 11.25.Wx, 14.40.Cs, 12.40.Yx 
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I. INTRODUCTION 

The nature of the light scalar mesons remains a key challenge for our understanding of 
how QCD generates bound states. Exceptionally strong and diverse mixing patterns between 
the differently flavored quarkonium components and with scalar glueballs are among the 
intricacies which kept this sector exciting and controversial for more than four decades. 
Theoretical approaches to the subject include a long tradition of model building, methods 
based on dispersion relations and effective field theories, QCD sum-rule analyses and more 
recently devoted and unquenched lattice simulations. The corresponding literature may be 
traced at least partly from the reviews lM5|. 



» 



In recent years, evidence has accumulated in favor of Jaffe's classic suggestion [6| that 
the lightest scalar meson nonet with masses below 1 GeV (and potentially a first recurrence 
bey ond 1 GeV [7]) may contain a dominant (cryptoexotic) tetraquark component [4J, |5|, [si 



1 if ] . This at present arguably most convincing interpretation demands an explanation for 
the exceptionally strong binding between the four "valence" quarks, however, which is able 
to reduce the lightest tetraquark mass significantly below that of the scalar quark-antiquark 
ground state. With the advent of holographic strong-coupling techniques based on the 
gauge/gravity correspondence 12[ and the__ construction of a first generation of approximate 



holographic QCD (or AdS/QCD) duals 13( it seems pertinent to explore whether those can 
provide new insights into existence and structure of strongly bound tetraquarks. In the 
present note we will take a first step in this direction. 



Several holographic studies of the light spin-0 meson [14J4l6l| (and glueball jl7|]) sec- 
tors in the AdS/QCD framework have recently appeared. Those were based exclusively on 
quark-antiquark interpolators, however, which presuppose a gg-dominated meson structure. 
Moreover, the first two of them exclude the lightest scalars (/o (600) and Kq (800)) from their 
analysis. In any case, it is rather obvious (and will be discussed in more detail below) that 
a straightforward extension of this approach to q 2 q 2 interpolators would describe four-quark 
states which are heavier, not lighter than the qq ground state. Our main focus will therefore 
be to search for a minimal but systematic AdS/QCD extension in which tetraquarks emerge 
naturally as the lightest scalar mesons. The sought-after holographic mechanism should 
further be able to render higher tetraquark excitations heavy (and thus broad) enough to 
avoid generating more low-lying scalar nonets than experimentally seen. 



In looking for a holographic representation of the required binding, we are guided by the 
expectation that the lightest scalar tetraquarks essentially consist of a "good" (i.e. maxi- 
mally attractive) diquark and a good antidiquark [38| which form an s-wave bound state 
pl, la [lQ|. The exceptional lightness of the tetraquark ground state is then mainly a con- 



sequence of the good diquark's large (and potentially instanton- induced |19l-l2l|) binding 



energy |22J. This suggests to adapt the approach introduced in Ref. [23|], which describes 
systematic mass reductions among the light baryons due to good-diquark binding holograph- 
ically in terms of bulk-mass corrections for their dual modes, to the problem. Below we will 
explore to what extent such corrections, as induced by the anomalous dimension of the 
corresponding four-quark interpolator, are able to generate light scalar tetraquarks as well. 

II. SCALAR SOFT-WALL DYNAMICS 

As a specific dynamical basis for our study we adopt the dilaton-soft-wall gravity dual of 



Ref. 24j which is governed by the action 

^ w) = ^/^V^^tr{|M| 2 -m^|X| 2 + i l f (F| + ^)}. (1) 

The gauge-covariant derivative Dm = du — ^Al,m + ^r,m contains the left/right fla- 
vor gauge potentials A L/R)M with field strengths (F L/R ) MN = d M A L/R)N - d N A L/RM - 
i [Al/r,m,A l /r,n]- The spin-0 bulk field 

X = (X + <p°t° + <p a t a ) e 2 "^ (2) 

consists of the background field X dual to the quark condensate, a flavor nonet ip A (with 
A G {0, ...8}) of scalar fluctuations (where the flavor singlet if may contain a glueball 
component) and the chiral pseudo-Goldstone bosons n. The U(3) flavor generators are 
normalized as tr \t A t B } = 5 AB /2. The remaining, non-dynamical background fields comprise 
the AdS 5 geometry 

R 2 
ds 2 = gMN<lx M ' dx N = — (rj^dx^dx" — dz 2 ) (3) 

with curvature radius R (where 77 is the four-dimensional Minkowski metric and z the 
Poincare coordinate of the fifth dimension) and the dilaton field 

$ (z) = X 2 z 2 (4) 



which is solely responsible for conformal-symmetry breaking in the infrared. The mass 
spectrum of the scalar meson nonets is then determined by the ^-dependent, bilinear part 

S[lf] = ^J SX >^ C ~* [9 MNd M^d N cp - m\^\ (5) 

of the bulk action ([T]). Since the dynamics (jSJ) is flavor diagonal, we have dropped the index 
of the scalar fields. (The "inverted" flavor-breaking pattern in tetraquark nonets [10| is thus 
beyond the scope of the present study.) 

Variation of the action fl5]) generates the field equation 

[V M V M - (<9 M $) d N + m%l<p = (6) 

where VmV 9 = ipM + ^mn) V 9 an d ^mn = ^mI* 1 \g\ ■ in terms of the four-dimensional 

Fourier transform 

d 4 q 

(^ 
the ensuing radial equation becomes 



if (x, z) = J -^je-^ip (q, z) (7) 



d 2 z -[- + 2\*z)d z + q 2 - R2ml 



Redefining the scalar field as 



z I z 2 



£^///2, 



<p(q,z) = 0. (8) 



<p{q,z)={-yer*-'*<l,(q,z) (9) 



then eliminates the first-derivative term and turns Eq. (JHJ) into the Sturm-Liouville problem 

[-d 2 z + V(z)]<j>{q,z) = q 2 cj>(q,z) (10) 

with the potential 

^) = (j + m^l + A 2 (AV + 2). (11) 



This formulation is convenient for evaluating the discrete meson mass spectrum m n = g„ 
associated with the normalizable scalar bulk modes <ft fan, z). 

At this point, a few comments are in order on the application of current AdS/QCD duals 
to three-color (i.e. N c = 3) physics in general and to exotics in particular. The semiclassical 
treatment of the underlying string (or, at weak curvature, gravity) dynamics is justified for 
large N c . Since calculating even the leading 1/N C corrections is beyond present capabilities, 
the extrapolation to N c = 3 has to remain naive and speculative. This holds in particular 



for bottom-up approaches where the N c dependence enters by comparison of calculated 



amplitudes with gauge-theory results 25 1 , for example in bulk coupling and mass parameters 
(one finds e.g. g\ = 12ir 2 /N c for the gauge coupling in Eq. (pQ)), in the normalization of 
correlators, or when encoding the Yang-Mills /3 function. From a formal perspective, the 
resulting extrapolation to smaller N c should therefore be regarded as a working hypothesis. 
In any case, it provides a consistent N c counting for the so far calculated observables and a 
semi-quantitative description of classical hadron properties (at the 10-30% level). 

One should keep in mind, however, that this naive treatment faces more challenging 
tasks in the exotic multiquark sector where the N c dependence can be particularly strong 
and where the extrapolation to N r = 3 becomes a practical necessity. Indeed, exotics are 



dissolving in the N c — > oo limit 26] (tetraquarks perhaps already for N c > 30 — 50, as 
unitarized chiral perturbation theory suggests [27|) and the higher- lying tetraquarks are 
expected to remain highly unstable under strong decay (with widths comparable to their 
masses) even at N c = 3 [6]. Our main focus will be on the lightest scalar tetraquark nonet, 
however, which should be relatively narrow (due to the paucity of open decay channels, 
limited phase space and partial dynamical suppression) and hopefully reasonably amenable 
to the naive N c — > 3 extrapolation. 

III. INTERPOLATOR AND ANOMALOUS DIMENSION 

The AdS/CFT dictionary [12j relates the normalizable bulk modes <p to those interpo- 
lating fields of the dual gauge theory which have enhanced overlap with the corresponding 
hadrons. More specifically, in the scalar sector the conformal dimension A of a meson 
interpolator prescribes the mass term in the potential (ITTj) as 

m\R 2 = A (A - 4) (12) 

and thereby sets the 2: — >• boundary condition for its dual mode, i.e. for the solution of 
Eq. ( fit)]) . In the case of ordinary scalar mesons the quark-antiquark interpolator 

4 = n A q a (13) 

is (up to factors) the unique choice with Aq q = 3. For light scalar mesons with a dominant 
tetraquark component, on the other hand, the larger number of possible couplings between 



the four quarks [28] allows for several potential tetraquark interpolators with the minimal 
classical scaling dimension A § 2 (? 2 = 6. In bottom-up duals it is therefore not a priori clear 
how to assign states to these interpolators. In the light of our above comments a natural 
choice, containing a good diquark and a good antidiquark, would be 

J* q2 = e abc e ade q b CT A q c q d CT A q e (14) 

where C is the charge conjugation matrix, no sum over A is implied and the s pin- flavor 
matrices T for the different members of the nonet may be found e.g. in Refs. [9|, |20|. The 
explicit form (fT4"]) is just an illustrative example, however, and will not affect our results. 
Indeed, for our purpose it suffices to specify the tetraquark interpolator's quantum numbers, 
its scaling dimension and the defining property of maximal overlap with the tetraquark 
ground state. 

The soft-wall meson spectrum resulting from the qq interpolator (1131) is [14j 



m^ n = 4(n+-)A 2 (15) 



(a) 



for n > 0. The eigenmodes are 4>q q (m n ,z) ex (Xz) {z/R)~ ' e~ x z ^ 2 L„ (A 2 z 2 ) where L 
is a generalized Laguerre polynomial 29] . The naive tetraquark spectrum corresponding to 
an interpolator like Eq. ( Tl4l) with classical scaling dimension A 5 2 g 2 = 6 (and no anomalous 
dimension), on the other hand, is 

m 2 A=6>n = 4 (n + 3) A 2 . (16) 

The radial excitations are thus likewise organized into a linear square-mass trajectory with 
the universal slope 4A 2 . The corresponding normalizable eigenmode solutions are now 

<M=6 K z) = iV A=6jn (A^) 6 (I) " 3/2 e~ x2z2 / 2 L^ (AV) . (17) 

The intercept of the A = 6 trajectory is twice as large as that of its A = 3 counterpart, 
however, which implies ^A = 6o/ m ggo = ^ fc> r the ground states. This exemplifies the general 
case: in the absence of bulk interactions and when ignoring anomalous dimensions (which is 
so far common practice in AdS/QCD applications), the holographic meson-mass predictions 
increase with the interpolator dimension A (as m 2 = 4 (n + A/2) A 2 in the dilaton soft-wall 
for d = 4). 



The above pattern adequately reproduces the empirical mass hierarchy among the less 
strongly bound qq mesons (with A = 3) and scalar glueballs (with A = 4). The exception- 
ally strong attraction required to reduce the masses of the lightest tetraquark nonet below 
those of the scalar qq ground state, on the other hand, is obviously not reflected in this 
minimal description. In the following we will examine whether the neglect of the tetraquark 
interpolator's anomalous dimension may be responsible for this shortcoming. As already 
mentioned, this would not be unexpected if the additional binding originates from the maxi- 
mally attractive diquark channel because diquark-content-dependent anomalous dimensions 
of the nucleon interpolators were found to describe good-diquark binding effects in the light 



baryon spectrum [23j (on the basis of the "metric soft- wall" gravity dual of Ref. |30|). In- 
deed, the anomalous dimension of the tetraquark interpolator will generally depend on its 
diquark content as well, and thus on the diquark content of the scalar meson with which 
it has enhanced overlap. Moreover, if two quarks form a good diquark whenever possible, 
and if good diquarks systematically reduce the masses of the hadrons which they dominate 
(compared to those containing other types of diquarks), one would expect this behavior to 
have a universal representation in the gravity dual. 

Holographically, the inclusion of the anomalous dimension 7 turns the mass term ffT2|) 
for the bulk modes dual to the A^2 g 2 = 6 interpolator into 

ml(z)R 2 = [6 + 1 (z)][2 + 1 (z)]. (18) 

(In more complex gravity duals, mass corrections may further arise from couplings to ad- 
ditional, e.g. condensate- related bulk fields.) The z dependence of 7 is inherited from its 
renormalization group (RG) flow 7 (/i) since the AdS/CFT dictionary translates the RG scale 
\i into the inverse of the fifth dimension, i.e. fi ~ 1/z [39]. Ideally, 7 (z) should be obtained 
from (both perturbative and nonperturbative) QCD and, as any other information from the 
gauge-theory side, implemented into the gravity dual in bottom-up fashion. However, at 
present this approach is ruled out by the absence of QCD information on 7 (z). As a prelim- 
inary substitute, we will therefore resort to a reasonable ansatz for 7 (z) and restrict it as 
tightly as possible by consistency, stability and physics requirements. This strategy cannot 
directly relate the ensuing mass corrections ( fl8l) to the anomalous dimension in QCD, but it 
will provide quantitative orientation and general insights, especially on how much tetraquark 
binding the anomalous- dimension mechanism can generate in principle, and even lead to a 



rather unique semi-quantative guess for 7 which is essentially independent of the specific 
AdS/QCD dynamics ([I]) adopted to derive it. Nevertheless, this estimate will have to be 
checked against - and in case of significant disagreements replaced by - QCD results for 7 
(e.g. from the lattice) when they eventually become available. 

A glance at the soft-wall potential (II ip with the generalized mass term (I18p for A 5 2 g 2 = 
6 reveals that even a scale-independent anomalous dimension 7 < — 3 would reduce its 
repulsion sufficiently to bring the tetraquark ground-state mass down below that of the qq 
ground state. However, it would do the same for all radial excitations as well, whereas 
phenomenologically higher-lying tetraquark excitations (with the exception of the first one 
if interpreted as a recurrence of the tetraquark ground state [7|) should be pushed up toward 
or beyond the corresponding qq excitations, to be broad enough to explain the absence of 
additional scalar nonets in the measured meson spectrum. This suggests to check whether 
a suitably running 7 (z) can accommodate the phenomeno logical situation. As it turns out, 
the necessary qualitative requirements may be cast into the simple power-law behavior 

7 (z) = -az v + bz K (19) 

with k ^> 7], so that the first (second) term modifies mainly the small-z (large-;?) region 
of the Sturm-Liouville potential. Moreover, the coefficients and exponents in Eq. (Tl9|) can 
be almost quantitatively determined by the underlying physics. To show this, we split the 
soft-wall potential ( TTTT) with the generalized mass term ( Tl8l) as 



V qH . (z) = (^ + 12^ i + A 2 (A V + 2) + AV (z) (20) 

into the uncorrected (i.e. 7 = 0) potential with the spectrum (fl6l) and the finite-7 correction 

AV(*) = 7 (*)[ 7 (z) + 8]i. (21) 

z z 

Maintaining asymptotically linear trajectories and enough repulsion to move the masses 
of the higher excitations upward restricts the exponent k to k < 2. Approximate slope 



universality 31] further demands |6| <C 1. The collapse of solutions to the mode equation 



(TTUj) into the AdSs boundary at z — ¥ 32J is avoided for 7] > (see below), and maximal 
attraction at smaller z, where the z~ 2 enhancement reduces the ground-state mass most 
effectively, demands 77 >C 1. As a consequence, Eq. (fl9l does not affect the AdS/CFT 
boundary condition of the tetraquark bulk modes for z — > since 7 (0) = 0, in accord with 
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asymptotic freedom. Maximizing the overall attraction of the potential correction (I2"T|) . 
finally, requires a ~ 4 (see below). It is worth noting that anomalous dimensions of a z 
dependence qualitatively similar to Eq. (I19p (in the region of interest) have recently been 
found in dual backgrounds which are of holographic RG-flow type [33] . (For a holographic 



estimate of the RG flow of an effective QCD coupling see Ref. |34j.) 

To put the above arguments in favor of an unusually large anomalous dimension for 
the tetraquark interpolator into perspective, we recall that in AdS/QCD applications to 
ordinary hadrons anomalous dimensions are generally neglected. This provides a decent ap- 
proximation as long as the anomalous dimensions of the corresponding two- and three-quark 
interpolators are much smaller than their classical scaling dimensions. In the UV, the latter 
is borne out by perturbative QCD results which show a weak, logarithmic RG scaling behav- 
ior. In the IR, on the other hand, where much less QCD information is currently available, 
it seems that anomalous dimensions can become relevant in some channels and generate 



notable holographic effects, as studied in Ref. [23]. Nevertheless, as a rule one does not ex- 
pect them to reach magnitudes comparable to the classical scaling dimensions. (Otherwise 
results of standard AdS/QCD would have to be revised). Our above arguments suggest that 
the exotic tetraquark channel provides an exception to this rule in which extraordinarily 
strong dynamical effects (related to good-diquark binding) render the anomalous dimension 
non-negligible even to leading approximation. 

IV. RESULTS AND DISCUSSION 

Several useful and rather generic features of the tetraquark mass spectrum in potentials 
of the type fl20|) . ( }2~Tj) can be understood without actually solving the eigenvalue problem 
( II Op . First, any 7 (z) which renders the correction (12 ip negative and decreases the power of z 
below —2 for z — > would cause a physically unacceptable "collapse to the center" instability 



of the Sturm-Liouville solutions 32]. In our context, the "center" at z = corresponds to 
the AdS 5 boundary, and this collapse is avoided by keeping 77 > 0. An instability may 
appear even if the singularity remains oc l/z 2 , however. Indeed, as long as the coefficient of 
this singularity is negative and exceeds a critical magnitude, the negative potential energy 
will overcome the stabilizing kinetic (or localization) energy. A straightforward analysis of 
the small- z behavior of the solutions of Eq. (llOp in the potential ( 12 Op for general boundary 



dimension d shows that such a collapse is prevented as long as 



- + A (A - d) > 0, (22) 



where we abbreviated A = A + 7. In fact, taking A constant turns the inequality (1221) into 



the Breitenlohner-Freedman (BF) bound m\R 2 > —d 2 /4 35]. A violation of this bound thus 
causes the ground-state wavefunction to disintegrate and its eigenvalue to become uiq = —00. 
This is a holographic signal for a vacuum instability in the boundary gauge theory. Our 
form ( l2Tf) of AV^, however, which arises from the anomalous-dimension correction ( fTHl) . has 
the lower bound 

AV(z)>-^ (23) 

z z 

(for d — 4 and all z). Hence it preserves the BF bound and prevents a collapse from the 
outset. The inequality ( 123]) is saturated for (constant) A = 2 or 7 = —4. This generates the 
lower bound 

m q 2 q 2 fl > m A=2 ,o = 2A (24) 

for the ground-state mass. Therefore 2A is the minimal tetraquark mass which can arise from 
the bulk mass term (TT8|) with any anomalous dimension. (A coupling of the tetraquark's 
dual mode to an additional bulk field may release the bound ( 1231) for z > and thus generate 
lower masses.) 

As a result, we have shown that a suitable ansatz for the tetraquark interpolator's anoma- 
lous dimension can indeed generate an extraordinarily large holographic binding energy 
which drives the lightest tetraquark mass from about 40% above down to around 20% below 
the mass m qq - t o = a/6 A (cf. Eq. ( 1151) ) of the qq ground state. To provide a concrete example 
for a 7 (z) of the form (fT9l) which results in a bulk-mode potential with nearly maximal at- 
traction for the ground state, we choose the coefficients inside the narrow range established 
by the criteria given above as a = 4, b = 0.05, 77 = 0.001 and k = 2. The resulting potential 
( 120]) with a typical soft-wall IR scale A = 0.425 GeV is plotted in Fig. [TJ (For comparison, 
we also show the potentials corresponding to the A = 3 and A = 6 interpolators without 
anomalous dimensions, as well as to the one with A = 2 which saturates the BF bound.) 
Numerical integration of Eq. ( TTOl) in the plotted potential indeed yields a ground-state mass 
close to the bound (1241) : m,g2 q 2 = 2.11A. Although the particular ansatz ( 1T9l) for 7(2) is 
not unique, Fig. [1] reveals that the bounds deduced above determine the behavior of the 
(almost) maximally attractive potential in a practically model-independent manner. 
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We now turn to the radial excitation spectrum of the tetraquark. The masses of the 
first five excited states (in the potential of Fig. [T|) are shown in Fig. [2j In this part of 
the spectrum, the square masses increase considerably faster than their asymptotic, linear 
trajectory 

ml n -±? 4nv / FTF ^ 4A 2 n (25) 

with almost universal slope 4A 2 . In particular, they start to surpass the corresponding 
masses on the qq trajectory around the second excitation level (and those on the uncor- 
rected A = 6 trajectory from about the fifth excitation onward). As intended, the higher- 
lying tetraquarks can therefore become broad enough to escape detection. Whereas a fully 
quantitative comparison of Fig. [2] with the experimental scalar mass spectrum as well as a 
detailed analysis of the numerically generated eigenmodes and decay constants would seem 
premature in our explorative framework, a few additional and probably rather robust fea- 
tures of the resulting spectrum are noteworthy. Interpreting the second scalar nonet around 
1.5 GeV as predominantly the lightest qq state would imply the phenomenological mass 
ratio rnq2 q 2 /mq q fi ~ 0.8/1.5 ~ 0.53 which is somewhat larger than what the tetraquark in- 
terpolator's anomalous dimension by itself can generate. (The alternative interpretation of 
the second nonet as a tetraquark recurrence [7[ , on the other hand, would raise the question 
of how to assign the lightest qq state.) The remaining part of the ground-state mass ratio 
may e.g. be supplied by the anomalous dimension of the qq interpolator or by couplings 
of the dual tetraquark modes to additional background fields. An adequate description of 
molecular-binding or meson-cloud effects could require to take bulk-mode interactions into 
account, moreover, and the resonance widths will further modify the masses. 

A reassuring feature of our anomalous-dimension-mediated tetraquark binding mecha- 
nism is that it will work in essentially the same fashion in any of the current AdS/QCD 
gravity duals. Indeed, 7 enters the bulk dynamics and hence the tetraquark mode poten- 
tial only through the mass term of Eq. (ITT]) which universally appears in all AdS/QCD 



duals 40J] and which is model-independently fixed by the AdS/CFT dictionary in Eq. ffT8|) . 
As a consequence, the anomalous- dimension- induced correction (121 p and its impact on the 
tetraquark spectrum are independent of the specific AdS/QCD dynamics under considera- 
tion. The 7 (z) determined above should therefore generate a light scalar tetraquark with 
substantially more massive radial excitations in any AdS/QCD dual. Nevertheless, it would 
still be interesting to explore the quantitative impact of the correction fT2T|) in other and 
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especially backreacted versions of AdS/QCD (such as Ref. [371] ). 

The holographic picture of light tetraquarks emerging from our above results can be 
extended and tested further, e.g. by implementing flavor-symmetry breaking and bulk-mode 
interactions, analyzing the ensuing decay channels and studying pertinent correlators (which 
may contain both qq and q 2 q 2 interpolators). It would also be interesting to investigate the 
behavior of exotic multiquarks in a topological expansion with large numbers of both colors 



(N c ) and flavors (Nf) where Nf/N c is kept fixed 36]. 



V. SUMMARY AND CONCLUSIONS 

We have discussed several foundational aspects of light scalar tetraquark physics from a 
holographic perspective. In particular, we have identified an AdS/QCD representation for 
the exceptionally large tetraquark binding energy in terms of resolution-dependent bulk- mass 
corrections for the corresponding dual modes. These corrections are argued to originate (at 
least in part) from the running anomalous dimension of the four-quark interpolator which 
encodes the diquark content of the tetraquark. As a consequence, they are independent of 
the specific AdS/QCD dynamics under consideration. The minimal achievable tetraquark 
mass is determined by the Breitenlohner-Freedman bound, furthermore, which prevents the 
tetraquark's dual mode from collapsing into the anti-de Sitter boundary. 

In order to provide a quantitative example for this binding mechanism, we have imple- 
mented it into the dilaton-soft-wall gravity dual for holographic QCD. The bulk mass cor- 
rections which generate nearly maximal attraction in the tetraquark channel are found to be 
largely determined by consistency and stability requirements, and the minimal tetraquark 
mass is fixed at twice the dilaton-induced infrared scale. As a result, the generated attrac- 
tion is indeed strong enough to render the tetraquark the lightest scalar meson, about 20% 
lighter than the qq ground state. To meet quantitative phenomenological expectations would 
probably require some additional contribution, however, which may originate e.g. from the 
anomalous dimension of the qq interpolator or from additional bulk fields and interactions. 
The higher tetraquark excitations turn out to be heavier than their qq counterparts, finally, 
so that their increasing widths (due to increasingly many open decay channels, a larger 
phase space and potential dynamical enhancements) should prevent them from generating 
more low-lying scalar resonances than experimentally seen. 
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FIG. 1: The Sturm-Liouville potential (J20J) of the tetraquark's dual mode for the anomalous 
dimension (|19l) with a = 4, b = 0.05, r/ = 0.001 and n = 2 (thick line). For comparison, we also 
plot the potential for the four-quark interpolator with A = 6, 7 = (dotted line) and for the qq 
interpolator with A = 3 (dashed line). In addition, we plot the potential for A = 2 (thin line) 
which saturates the Breitenlohner-Freedman bound. 
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FIG. 2: The ground state (n = 0) and the first five excitations of the tetraquark mass spectrum in 
the potential of Fig. CD (thick bars). We further plot the spectra for the qq interpolator (Eq. (fT5J) . 
dashed bars) and for the A = 6 interpolator (Eq. (|16p . dotted bars) with vanishing anomalous 
dimensions. 
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